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In this paper, I relate the phase separation transition in binary hard-core mixtures in the limit of 
small size ratio to a bond percolation transition. It allows to estimate the critical packing fraction 
at the transition as a function of the size ratio and the composition of the mixture for different 
shapes of objects. The theoretical predictions are in excellent quantitative agreement with numerical 
simulations of binary parallel hard squares mixtures. 
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Binary mixtures of impenetrable objects pose an eas- 
ily formulated problem which is of crucial importance to 
understand entropic effects in colloidal systems. In this 
systems of "large" and "small" objects, an attractive de- 
pletion force between the large particles results from the 
presence of the small ones jl] . The main interest in such 
mixtures is the existence of a phase separation induced 
by this effective attraction of large objects. 

Recent Monte-Carlo simulations [^]]|] have shown un- 
equivocally that the phase separation takes place for 
spheres, parallel cubes and parallel squares for sufficiently 
dissimilar objects. Experiments Q| also predict phase 
separation for mixtures of hard spheres. Now that esti- 
mates for the critical packing fraction at the transition 
have become available, theoretical work has become all 
the more interesting. Closures of the Ornstein-Zernickc 
integral equations have been studied extensively, and nu- 
merical calculations with the Rogers- Young approxima- 
tion predict a demixing transition for sufficiently small 
size ratios ||. The main drawback of this approach is 
that the closures are essentially uncontrolled. Other ap- 
proaches have been used ||f7| which lead to an instabil- 
ity between a homogeneous fluid phase and a demixed 
fluid-fluid phase. However, as shown by numerical sim- 
ulations [g,D and experiments [Q], for small size ratio, 
the transition is between a fluid phase and a fluid-solid 
phase where the solid phase is constituted by closely 
packed large objects surrounded by a fluid of small ones. 
The actual nature of this "solid" phase is yet hotly de- 
bated ff-§. The aim of this paper is not to give a better 
understanding of this phase but to determine the critical 
packing fraction at the transition. 

In this article, I propose an expansion for the large- 
large pair distribution function in the fluid phase which 
is expected to be convergent in the limit of small size 



ratio. If two large particles are sufficiently close so that 
no room is left to insert a small particles between them, 
they are called bonded. This expansion, a virial develop- 
ment, allows to calculate the number of bonds between 
large-large objects, which is an increasing function of the 
packing fraction. Then, I relate the phase separation 
transition to a bond percolation transition and the solid 
phase to the percolating aggregate. Even if this heuristi- 
cal mapping does not give a better understanding of the 
solid phase, it allows to compare the theoretical predic- 
tions of the critical packing fraction and its dependence 
on the composition with numerical simulations |^,|3| and 
experiments |I| for objects of different shapes. 

In the following, I consider large and small spheres 
(or disks) of diameters o~i and o~ s and parallel cubes (or 
squares) of sides 07 and o~ s . The total packing fraction 
is i] = 7] s + r\i with r\ a = N a v a /V, a = s,l where V is 
the box volume, N a the number of particles and v a the 
volume of a particle. The size ratio and the composition 
are defined as R = o~ s /o~i and x = T\ijr\ respectively. For 
simplicity, distances are considered in units of the large 
particle size (07 = 1 so that a s = R). 

The pair distribution functions are the relevant func- 
tions to characterize the homogeneous fluid phase. The 
Mayer functions f a b allow to obtain density expansions of 
the large-large pair distribution function gu Jl(],[ll| . For 
cubes (d — 3) and squares (d = 2): 

d 

fab(r) = -l[e(a a + a b -2\n\) (1) 

i=l 

where rj stands for the difference of the i th component of 
the positions of particles a and b and 8 is the Heaviside 
step function. For spheres and disks: 

/ o6 (r) = -Q(a a + a b - 2r) (2) 
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where r is the distance between the particles a and b. The 
functions f a b are simply step functions with /„& = — 1 if 
the two particles a and b overlap and = if not. 

I propose the following exponential expansion of the 
pair distribution function which is a direct generalization 
of the one-component expansion (see Chap. 5 in 

gu(r) = exp [ r, s R- d a - 1 J\+ m^ 1 /\ 

* a < 2 n +4 n +h ) + •••] 

where d is the dimension of the space and a is a constant 
that depends on the shape of the objects (a = 7r/6 for 
spheres, a — ir/4 for disks and a = 1 for cubes and 
squares). Each diagram is an integral over positions 
of small and large particles represented respectively by 
small and large black dots. Each link corresponds to a 
Mayer function and the two large open circles to two 
large particles with position difference r. For example, 
the first diagram corresponds to: 

A = ffsi(r')fsi(r-r')dr> 

and is a simple function of the position r and the size 
ratio R. The diagrams in the expansion are all distinct 
connected diagrams consisting of two non-adjacent open 
circles, at least one black dot and such that the two open 
circles do not form an articulation p"T[ . A diagram with 
articulations may be reduced to a product of simpler ones 
(see below). 

The above expansion, a density development, has neg- 
ative powers of R when expressed with the pertinent pa- 
rameters: R, x and rj. However, three different effects 
may compensate the negative powers: first of all, a link 
between two small particles scales as R d since the cor- 
responding Mayer function restricts the integration to 
a small particle volume; second, diagrams of the same 
order but with odd and even numbers of links may com- 
pensate each other because the Mayer functions are neg- 
ative; and last one, the critical packing fraction at the 
transition scales as a positive power of R. Thus, we may 
expect that the expansion proposed in this article is con- 
vergent when written in powers of the size ratio at the 
phase separation transition. In that case, only the first 
few diagrams have to be considered in the limit of small 
R. 

In the article of Dijkstra et al. another expansion 
up to the second order was presented for gu . This expan- 
sion is not convergent in the limit of small R due to the 
presence of diagrams which form an articulation, like the 
one with two small particles without link between them. 
In fact, this diagram is the square of a simpler one: 



M = ( A) 2 - 



It scales with a smaller power of R than the others of the 
same order and corresponds to the most divergent part 
in the second order expansion (see Eq.(5) in [[10|). 




FIG. 1. This figure represents the excluded volume for 
large disks around the disk in gray. Large disks are bonded to 
the disk if their centers stay in the ring of width equals to the 
small particle diameter around the excluded volume. In this 
case, the disks 1 and 2 are bonded to the disk whereas the 
disk 3 is not. This case corresponds to a size ratio R — 1/20. 

In the binary mixtures, gu shows a sharp peak at the 
contact value. In the limit of small R, this peak is as 
narrow as R and corresponds to large particles with a 
distance 1 < r < 1 + R between them. Two such large 
particles may be considered as bonded. In this case, the 
mean number n b of bonds per large particle is given by 
the integration of gu on the shell of width R around the 
excluded volume multiplied by the density of large par- 
ticles pi = Ni/V: 

n b = pi / 9u(r)dr (3) 

Jl<r<l+R 

where r is the usual norm of r for spheres and disks, and 
it is the maximum of the absolute values of the coordi- 
nates of r for cubes and squares. In the latter case, it is 
the relevant distance since I consider only parallel cubes 
(or squares). The example of disks is represented in Fig.l 
where the gray region stands for the excluded volume of 
large disks around the disk 0. The ring around the gray 
region corresponds to the shell on which the integration 
is made. Centers of disks 1 and 2 are in this region and 
those disks are bonded to the disk but the disk 3 is not. 
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From gn , it is clear that the mean number rib of bonds 
per particle is an increasing function of the packing frac- 
tion. Thus, on increasing rj, larger and larger aggregates 
of large particles exist in the system. There exists a crit- 
ical number of bonds per particle n c for which a perco- 
lating aggregate appears. This corresponds to the phase 
separation transition and the percolating aggregate con- 
stitutes the solid phase. Notice that above this phase 
separation transition, gu is no longer defined since the 
translation invariance disappears due to the formation of 
the solid phase. If we suppose that bonds are randomly 
distributed, the critical number n c may be related to the 
bond percolation threshold p c of the lattice corresponding 
to the close packed configuration. The latter corresponds 
to the fraction of bonds necessary to obtain a percolating 
aggregate whereas n b counts the mean number of bonds 
per particle. Since the number of bonds per particle in 
the close packed configuration is the coordination num- 
ber z, n c = zp c . The fact that positions of objects are not 
restricted to lattice sites may slightly modify the critical 
number n c . However, z is an obvious upper bound for 
the critical number n c . 

The positions of close packed spheres describe a FCC 
lattice of coordination number z = 12 and the bond per- 
colation threshold of this lattice is p c — 0.199 p2| . Thus, 
the critical number of bonds per particle is n c = 2.4. 
Disks are close packed in a triangular lattice (p c = 0.347 
and z = 6) II which leads to n c = 2.1. Cubes and 
squares in a cubic or square lattice (p c = 0.249 or 0.5 
and z = 6 or 4 respectively) Jl^] have n c = 1.5 or 2 
respectively. 

In order to obtain the behavior of the critical packing 
fraction in the limit of small R, I consider only the first 
diagram of the exponential expansion assuming that the 
others are negligible. In this case, n& is given by: 

n b ^Jdv exp J dr' f sl (v')f sl (r - r')) (4) 

where the integration over r is restricted to 1 < r < 1+R. 
The critical packing fraction rj c is defined by rib{r} c ) = n c 
and depends on the size ratio and the composition. An 
upper bound rj up for this critical packing fraction is given 
by the condition nb(ji up ) = z. The estimate for r\ c is 
valid only for small R, for which the peak at the con- 
tact value of the pair distribution function of large-large 
particles may be related to bonds between large particles 
and where the diagram considered gives the dominant 
part of rib- 

Consider first parallel cubes and squares. The integra- 
tion of Mayer functions f s i is: 

„ d 

/ dv'.f sl {v').f sl {v-v') = IJ(l+ J R-|r i |)e(l+ J R-|r li |). 

(5) 



For cubes, the mean number rib deduced from this 
equation and Eq.(Q) leads to the following critical pack- 
ing fraction: 

, c (Cubes)_— log^j. (6) 

The upper bound for the critical packing fraction is given 
by the same equation replacing n c by z. Notice that n c 
used to calculate r\ c is not essential to obtain the asymp- 
totic behavior — 3i? 2 (l — x)^ 1 logi? which is also given 
by the upper bound r] up . The packing fraction of small 
cubes rj s = r}(l — x) at the transition is slightly dependent 
on composition (through a logarithmic term negligible in 
the limit of R — > 0). This result, which is observed for 
objects of all shapes, is not surprising since, in a first ap- 
proximation, only the small particles are responsible of 
the effective attraction of the large ones. 

The same calculation for the critical packing fraction 
may be done for parallel squares: 

l^^ 5 ^)' CD 

It is interesting to notice that this is the first theoretical 
prediction of such a phase separation transition in binary 
mixtures of parallel hard squares. 

The behavior of rj c as a function of R and x has been 
obtained for the first time in recent simulations of mix- 
tures of parallel hard squares ||. Thus, it is possible to 
compare the theoretical predictions with the direct sim- 
ulations. For this case of parallel hard squares, I have 
been able to push the expansion considerably farther (to 
second order) since the diagrams are quite simple to cal- 
culate. Mayer functions are products of one-dimensional 
functions and the same holds true for the diagrams. In 
the following table, theoretical predictions with first or- 
der and second order expansion are compared with sim- 
ulations. 



R 


X 


1 st order 


Vc 

2 nd order 


Simulation 


1/10 


0.3 


0.68 


0.52(0) 


0.49 ±0.02 


1/10 


0.5 


0.78 


0.54(4) 


0.53 ±0.02 


1/10 


0.7 


> 1. 


0.58(8) 


0.60 ±0.05 


1/20 


0.5 


0.53 


0.42(4) 


0.41 ±0.02 


1/50 


0.5 


0.29 


0.26(4) 


0.26 ±0.02 



For the first order, a qualitative agreement is obtained 
for R = 1/10 and 1/20, whereas for R = 1/50, the agree- 
ment is even quantitative. The agreement between nu- 
merical simulations and predictions with the second order 
expansion is quite striking, and indicates that this expo- 
nential expansion succeeds in summing up the essential 
diagrams. It is interesting to notice that, replacing the 
packing fraction by its critical value (Eq. (|7j) ) in the expo- 
nential expansion, the second order terms scale at least 
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as R whereas the first order term is of order one (ne- 
glecting logarithmic terms). It would be interesting to 
prove that higher orders are also negligible. Moreover, 
the width of the shell on which the integration is made 
to calculate n& does not seem to be essential as soon as we 
consider a width larger than the peak of the pair distri- 
bution function (~ R) but still of order R. I have checked 
for R = 1/50 that calculations with shells of width 2R 
and 3R do not modify significantly ij c . 

Let me now consider the mixture of spheres. The in- 
tegration of the Mayer functions has already been calcu- 
lated 0: 

J dr'f sl (r')f sl (r-r') = ^(l + R-r) 2 (l + R+r/2) (8) 

for 1 < r < 1 + R. Including this expression in Eq.ffl), I 
deduce the critical packing fraction: 



7] c (Spheres) 



2R 



3(l-x) 



log 



3 (1 — x) n c 
8ttxR 2 



(9) 



which is valid for r\ c 3> 2 i?/(3 (1 — xj). This implies that 
the logarithm of the size ratio is sufficiently large since 
T] c decreases to zero with R not only linearly but with 
logarithmic corrections. For x = 1/2 and R — 1/30, 
Eq.(Q) gives r\ c — 0.24(8) whereas numerical simula- 
tions have shown that r\ c < 0.243. Furthermore, both 
predictions for the critical packing fraction agree quali- 
tatively with experimental results Unfortunately, in 
order to be able to compare quantitatively theoretical 
with experimental predictions, either new experiments 
with smaller size ratios or theoretical predictions consid- 
ering further diagrams in the expansion beyond the first 
one are needed. 

An analogous calculation may be done for disks in two 
dimensions, where the first diagram scales as i? 3 / 2 in the 
limit i?— >0forl<r<l + i?. This leads to a critical 
packing fraction: 



i] c (Disks) 



3tt VR 
8V2 1 -x 



log 



'V2{l-x) nc 
2ir x R 3 / 2 



(10) 



The scaling prediction in V R for r\ c has been already sug- 
gested in H , but there is not yet numerical evidence for 
the phase separation transition. From Eq. (|i"o|) , in order 
to have r\ c ~ 0.6 at a composition x = 0.5, the size ratio 
to consider is R ~ 1/600. This seems yet beyond reach 
since for this size ratio, we would have to take into ac- 
count 360 000 small disks for each large one or, in order 
to clearly see the transition, a few dozen of large parti- 
cles and more than 10 7 small ones. On the other hand, 
taking into account higher orders in the expansion, the 
estimation of the size ratio may be strongly modified. 

In conclusion, an exponential expansion of the pair dis- 
tribution function of large-large particles in binary mix- 
tures has been presented. A known limitation of such 



expansion is of course that it is limited to the homoge- 
neous phase. In this work I have shown that, neverthe- 
less, this expansion can be used to predict a fluid-solid 
phase separation by computing the quantity which trig- 
gers it: a large increase in the number of large particles 
in a thin shell around a given large particle. This was 
already pointed out in [0J3]. In the limit of small size 
ratio, the mean number of bonds between large parti- 
cles at the phase separation transition is related to the 
bond-percolation threshold of the lattice corresponding 
to the close packed configuration, and the solid phase 
corresponds to the percolating aggregate of large parti- 
cles. The critical packing fraction is then deduced. Its 
dependence with the composition confirms that the effec- 
tive attraction of large objects is principally due to the 
depletion of small ones. Additional numerical simulations 
of mixtures of parallel hard squares reported elsewhere || 
confirm the theoretical predictions. The theoretical re- 
sults for mixtures of spheres are in qualitative agreement 
with experiments, whereas the phase separation transi- 
tion for mixtures of hard disks is predicted to occur at 
such small size ratio that the corresponding simulation 
above the transition seems yet out of reach. 

It is a pleasure to thank W. Krauth for illuminating 
discussions and M. B. Gordon and W. Krauth for a care- 
ful reading of the manuscript. 
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